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I971>  was  prepared  by  the  Westinghouse  Electric  Corporation.  Special 
recognition  is  due  Dr.  R.  C.  Smith,  of  the  Research  Laboratories,  who 
worked  out  the  theoretical  analysis  in  Section  4. 

The  research  reported  in  this  document  is  supported  by  the  Advanced 
Research  Projects  Agency  of  the  Department  of  Defense  and  is  monitored  by 
the  Bureau  of  Mines  under  Contract  HO1 10377. 
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SUMMARY 


PURPOSE 

Westinghouse  is  seeking  field  operating  data  with  which  to  verify  and 
extend  theoretical  and  laboratory  data  on  the  feasibility  of  employing  the 
electron  beam  gun  (EBG)  to  satisfy  the  needs  for  removal  of  hard  rock. 
These  practical  data  will  provide  initial  definition  of  equipment  con¬ 
figuration  and  applicable  modes  of  operation. 

TECHNICAL  PROBLEMS 

Support  Equipment.  Westinghouse  is  modifying  its  existing  36-kW  EBG  and 
will  provide  support  equipment  necessary  for  field  operation. 

Field  Tests .  This  program  involves  the  use  of  the  36-kW  EBG  to  determine 
whether  laboratory  test  data  on  melting  rates,  breakage  mechanisms,  etc., 
can  be  extrapolated  to  predict  field  performance.  The  field  site  being 
selected  will  provide  working  faces  representative  of  a  large,  relatively 
unfractured  rock  mass.  During  field  tests,  Wes‘:i  nghouse  wi  1 1  determine 
the  techniques  and  procedures,  such  as  cutting  pattern,  gun  tracing  speed, 
and  method  of  cut,  required  to  produce  the  optimum  excavation  rate. 

Theoretical  Studies.  Westinghouse  is  studying  the  effects  of  different 
heating  geometries  and  heating  rates  on  stress  distribution  patterns  and 
associated  breakage  phenomena.  Westinghouse  will  also  attempt  to  develop 
theoretical  models  for  an  optimum  mode  of  operation. 

Laboratory  Tests.  Westinghouse  is  conducting  laboratory  tests  to  support 
theoretical  work  and  field  experiements  and  to  extend  the  use  of  the 
EBG  to  different  rock  types. 

Systems  Analysis  Study.  Westinghouse  will  evaluate  the  feasibility  of 
developing  a  hard  rock  excavation  system  using  the  EBG  as  the  prime  source 
of  power  for  rock  fragmentation. 
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TECHNICAL  RESULTS 

This  report  covers  effort  performed  during  the  second  quarter  of  the  con¬ 
tract.  This  work  is  a  continuation  of  that  reported  i.n  the  first  quarterly 
report  (Ref.  3)’f  Work  was  performed  in  four  of  the  five  areas  summarized 
abovi;  the  systems  analysis  study  will  not  become  active  until  the  field 
test  effort  has  been  completed  and  evaluated.  A  sixth  area  (not  included 
in  the  contract)  is  the  completion  of  the  36-kW  EBG  and  its  conversion 
from  a  laboratory  model  to  a  field  test  version. 

Support  equipment.  Design  work  and  manufacturing  effort  for  the  EBG 
carriage  were  completed  in  June  1971*  The  carriage  is  presently  under¬ 
going  a  series  of  preliminary  operational  and  calibration  tests.  Some 
deficiencies  in  carriage  design  were  noted  in  these  tests,  but  these  are 
being  corrected  prior  to  installation  of  the  36-kW  EBG.  Additional  tests 
will  be  conducted  with  the  EBG  installed  and  with  operational  field  support 
equipment  to  verify  system  operation. 

Field  Tests.  Field  test  activity  in  the  second  quarter  of  the  contract  has 
again  been  primarily  in  the  area  of  site  review  and  selection.  The  site 
previously  selected,  i.e.,  the  Logan  quartz-gabbro  quarry,  is  no  longer 
being  considered.  The  Logan  rock  is  extensively  fractured  and  faulted  and, 
therefore,  was  deemed  geologically  undesirable,  although  it  met  all  other 
criteria  for  site  selection.  Another  site  has  been  selected,  which  fulfills 
the  geological  requirements  and  most  other  site  criteria.  This  is  the 
Raymond  Granite  Company's  quarry  in  Raymond,  California,  which  is  about 
25  miles  south-southwest  of  Yosemi te  National  Park. 

Other  activities  in  the  area  of  field  tests  have  been  the  preparation  of  a 
preliminary  field  test  plan,  the  specification  for  and  concept  design  of  a 
radiation  shield  for  field  use,  and  specification  and  acquisition  of  other 
site  support  equipment,  such  as  motor-generator  systems,  radiation  moni¬ 
toring  and  detection  instrumentation,  and  the  required  site  utilities. 

*See  list  of  references  in  Section  6. 
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The  EBG  as  modified  for  rock  cutting  In  the  field  (long-beam  transfer 
section  to  reach  into  holes  and  crevices)  is  operational  at  low  power. 

Some  delays  were  encountered  in  going  to  the  desired  higher  power  levels, 
and  consequently  shipment  to  the  field  site  was  delayed.  This  should  be 
accomplished  in  k  to  6  weeks. 

Theoretical  Studies.  A  mathematical  model  has  been  developed  that  can  de¬ 
scribe  the  piercing  of  the  rock  by  the  electron  beam  and  the  thermal  stresses 
which  are  subsequently  produced. 

Computer  printouts  have  been  obtained  for  the  stresses  after  5  to  300 
seconds  of  piercing  with  a  9~kW  beam.  It  is  not  easy  to  derive  from  these 
plots  when  and  where  the  rock  will  crack,  but  trends  in  the  time  develop¬ 
ment  of  the  stress  field  can  be  seen  that  confirm  laboratory  observations, 
at  least  for  one  particular  geometry.  Computations  must  now  be  extended 
to  other  types  of  rocks,  power  levels,  and  geometric  conditions. 

Laboratory  Tests.  No  laboratory  tests  were  performed  during  the  second 
quarter. 

DoD  IMPLICATIONS 

The  primary  DoO  implication  of  a  successful  field  test  program  is  the 
potential  to  fragment  and  remove  hard  rock  from  any  location  without  de¬ 
grading  the  structural  integrity  of  the  parent  rock  from  which  it  was 
removed. 

Upon  completion  of  the  appropriate  development  program,  the  EBG  system 
could  be  applied  to  rock-removal  systems  for  programs  such  as 

Sangui ne 

Cheyenne  Mountain 
Mi nuteman 
Safeguard 
AEC  testing 

•  • 
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and  probably  a  number  of  other  applications.  In  addition,  an  operational 
EBG  system  would  have  applications  to  the  programs  of  the  Department  of 
Transportation  and  other  agencies. 

It  should  be  stressed,  however,  that  a  considerable  development  program 
will  be  required  after  the  assumed  successful  field  tests  before  an 
operational  system  can  be  delivered. 

IMPLICATIONS  FOR  FURTHER  RESEARCH  AND  DEVELOPMENT 

Westlnghouse  believes  that  the  current  field  test  program  will  demonstrate 
the  potential  of  the  EBG  for  hard-rock  fragmentation.  Based  on  the  success¬ 
ful  completion  of  these  tests,  additional  research  and  development  effort 
can  be  envisioned  In  the  following  areas: 

1.  Additional  tests  In  various  types  of  open  quarrries  and  in 
underground  mines 

2.  Evaluation  of  EBG  in  conjunction  with  existing  mining  and  tunneling 
equipment  and  as  part  of  a  total  system  design. 

3.  Improved  knowledge  of  rock  failure  mechanisms. 

k.  Investigation  of  the  use  of  ion  beams  with  their  reduced  shielding 
requirements. 

5.  Increase  in  power  and/or  Increase  In  electron  acceleration 
vol tage. 

6.  Real-time  predictions  of  geologic  conditions  using  existing  optical 
and  acoustic  sensing  technologies. 

These  areas  of  continuing  investigation  will  augment  a  product  development 
program  that  is  aimed  at  the  achievement  of  a  system  to  satisfy  govern¬ 
ment  and  industrial  needs  for  hard-rock  removal. 
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1.  INTRODUCTION 

Laboratory  tests  conducted  by  Westinghouse  have  shown  that  the  electron 
beam  gun  (EBG)  may  be  a  promising  tool  for  use  in  hard-rock  excavation 
(Refs.  1  and  2)*.  These  tests  have  been  conducted  on  small,  unconstrained, 
laboratory  test  specimens  weighing  a  few  hundred  pounds.  Under  the  present 
contract,  field  tests  are  to  be  conducted  to  study  the  electron  beam 
process  on  a  semi -i nf i ni te,  constrained  rock  mass,  such  as  might  be  en¬ 
countered  in  mining  or  tunneling  operations. 

The  objectives  of  the  electron  beam  gun  evaluation  program  are: 

1.  To  obtain  field  operating  data. 

2.  To  determine  the  effectiveness  and  economic  feasibility  of 

the  electron  beam  gun  compared  with  conventional  methods  of  hard- 
rock  excavation. 

3.  To  determine  practical  optimums  for  equipment  configuration  and 
modes  of  operation. 

In  support  of  the  foregoing  effort,  laboratory  experiments  and  theoretical 
studies  are  being  conducted  to  compute  and,  if  possible, to  predict  the 
thermal  stresses  and  the  resultant  rock  fragmentation  for  various  cutting 
strategies  and  electron  beam  parameters  in  different  types  of  rock. 


*See  list  of  references  in  Section  6. 
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2.  ELECTRON  BEAM  GUN  AND  SUPPORT  EQUIPMENT 

2.1  ELECTRON  GUN  AND  CONTROLS 

The  electron  beam  gun  ( EBG )  was  described  in  the  First  Quarterly  Report 
(Ref.  3).  At  this  time  the  gun  has  been  operated  at  150  kv  and  a  few 
kilowatts  of  beam  power,  but  adjustment  of  the  electron  optics  for 
beam  guidance  through  the  long-beam  transfer  section  at  high  power  has  not 
been  completed.  Unfortunately,  this  has  delayed  the  start  of  the  field 
tests . 

2.2  SUPPORT  EQUIPMENT 

Supprrt  equipment  for  the  EBG  consists  of  the  gun  carriage,  the  vacuum 
system,  and  the  electrical  power  supply  system.  Except  for  the  motor- 
generator  sets  required  in  field  operations,  the  vacuum  system  and  elec¬ 
trical  system  existed  as  laboratory  equipment  prior  to  start  of  work  under 
this  contract.  Thus,  the  principal  task  in  this  area  is  the  design  and 
manufacture  of  a  carriage  to  support,  position,  and  manipulate  the  gun 
during  field  test  operations.  The  design  task  was  completed  in  April  1971. 
A  complete  design  description  appeared  in  the  First  Quarterly  Report 
(Ref.  3,  p.  2-6).  The  manufacturing  portion  of  the  task  was  begun,  and 
the  completed  carriage  was  delivered  to  Westinghouse  in  mid- June  1971. 

Work  on  the  carriage  since  delivery  has  consisted  of  completing  the  elec¬ 
trical  and  control  systems,  conducting  an  operational  test  of  the  carriage, 
and  making  minor  modifications  as  a  result  of  deficiencies  observed  in  the 
tests.  These  tests  and  modifications  are  described  below. 

2.2.1  Ca rriaqe  Tests 

In-plant  testing  of  the  EBG  carriage  Is  subdivided  into  two  categories: 

(1)  an  initial  phase  using  a  mock-up  of  an  EBG,  and  (2)  a  final  phase  to 
verify  operation  of  the  carriage  in  conjunction  with  the  36-kW  EBG. 
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The  initial  phase  of  the  in-plant  tests  consists  of  several  tests  to  verify 
operation  of  the  carriage  hydraulic  system,  to  obtain  calibration  data  for 
the  system  velocity  controls,  and  to  detect  and  measure  any  beam  oscilla¬ 
tions  that  might  occur  as  a  result  of  hydraulic  transients  or  structural 
vibrations.  Preliminary  carriage  tests  were  performed  in  a  no-ioad  condition 
to  verify  operation  of  ail  components.  A  concrete  and  steel  pipe  mock-up 
of  the  EBG  was  then  constructed  and  mounted  on  the  carriage  for  the  remain¬ 
ing  tests.  The  mock-up  gun  approximates  the  field  test  gun  in  space  envelope, 
weight,  and  center-of-grav i ty  location.  The  carriage,  with  the  mock-up  gun 
installed,  appears  in  Figures  2-i  through  2-3.  Figure  2-3  is  a  time  exposure 
taken  during  carriage  translation  with  a  light  attached  to  the  end  of  the 
transfer  column  to  record  the  gun  motion.  Although  ail  modes  of  operation 
are  represented  in  Figure  2-3,  these  motions  do  not  extend  to  the  limits 
of  the  gun  manipulator.  The  control  pane!  for  the  carriage  is  shown  in 
Figure  2-4i . 

After  installation  of  the  mock-up,  a  series  of  calibration  tests  was  run 
to  obtain  the  data  required  for  velocity  control.  Each  portion  of  the 
carriage  hydraulic  system  is  equipped  with  temperature-  and  pressure- 
compensated  variable  flow  control  vaives.  These  vaives  wiii  permit 
repeatable  and  constant  velocity  control  for  each  mode  of  operation.  The 
calibration  consists  of  operating  the  various  portions  of  the  hydraulic 
circuit  at  several  flow  control  settings.  The  actuation  time  for  each 
operation  is  recorded  and  the  results  plotted  as  beam  velocity  vs.  valve 
setting.  These  tests  will  be  completed  In  eariy  Juiy. 

The  second  phase  of  testing  wiii  begin  when  the  36-kW  gun  is  installed  on 
the  carriage.  This  phase  is  comprised  of  a  series  of  tests  to  recheck  the 
previous  valve  calibration  data,  to  identify  any  physical  interference 
problems  that  may  exist,  and  to  verify  operation  of  the  EBG  with  the  field 
configuration  of  vacuum  pumps  and  motor-generator  sets.  The  operational 
tests  wiii  include  some  preliminary  cutting  tests,  to  be  performed  using 
a  6  x  6  x  5-foot  granite  test  block. 
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F igure  2-1.  Carriage  with  Hock-up  of  Electron  Bear,  Cun  (Front  VI 


Figure  2-2.  Carriage  with  Mock-up  of  Electron  Beam  Gun  (Right  Side  View) 
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2.2.2  Mod i f icat ions 

In  the  preliminary  testing,  some  discrepancies  in  carriage  operation  were 
noted.  The  clearances  between  the  portions  of  the  gun  mount,  which  move 
in  the  Y-axis  or  vertical  mode,  were  found  to  be  excessive,  and  under 
certain  loading  conditions,  the  movable  portions  could  be  "derailed."  A 
system  of  auxiliary  guide  rollers  has  been  designed  and  installed  to  redis¬ 
tribute  the  loads,  and  thereby  alleviate  the  problem.  As  mentioned  in  the 
First  Quarterly  Report  (Ref.  3),  a  flexible  coupi ing  was  provided  in  the 
gear  drive  system  to  prevent  any  misalignment  problems  that  might  result 
from  a  redundant  center  condition.  Prior  to  its  installation,  the  gear 
drive  coupling  was  a  rigid  connection.  The  rigid  coupling  appeared  to 
function  properly  but  was  removed  anyway  and  replaced  with  the  flexible 
device.  The  flexible  rubber  insert  in  the  new  coupling  resulted  in  soft 
or  spongy  response  in  the  yaw  drive  with  a  considerable  amount  of  over¬ 
shoot  observed  at  the  end  of  the  actuator  stroke  and  upon  controi  valve 
closure.  The  rigid  coupling  will  be  reinstalled  to  eliminate  these  un¬ 
desirable  characteristics. 
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3.  FIELD  TESTS 

3.1  TEST  SITE  SELECTION 

Site  review  and  selection  have  been  conducted  with  considerations  such  as 
site  geology,  availability  of  support  equipment  and  facilities,  and  location 
in  mind.  Toward  the  end  of  the  first  quarter  of  this  contract,  a  tentative 
site  selection  was  made  after  reviewing  12  quarry  sites  and  an  underground 
mine.  The  site  selected  was  the  Logan  quartz  gabbro  quarry,  a  producer  of 
crushed  rock  products,  located  near  Watsonville,  California.  This  site  met 
most  of  the  criteria  for  an  acceptable  site  except  the  requirement  for 
relatively  unfractured  rock.  Since  it  is  located  on  the  San  Andreas  Fault, 
the  Logan  rock  is  extensively  fractured  and  :aulted.  After  a  second  visit 
to  the  Logan  site  by  Westinghouse  and  Bureau  of  Mines  personnel  in  June  1971, 
this  site  was  rejected  as  a  prime  test  site,  and  another  site  was  located 
that  more  completely  fulfills  the  desired  geologic  requirements. 

The  approved  test  site,  also  visited  by  Westinghouse  and  Bureau  of  Mines 
personnel,  is  the  Raymond  Granite  Company's  "Sierra  White"”  granite  quarry. 
The  quarry  is  located  In  the  Sierra  Nevada  foothills  of  California  and 
25  miles  south-southwest  of  Yosemite  National  Park.  A  section  of  the 
USGS  map  for  California  gives  an  indication  of  geological  formations  In 
the  Raymond  area.** 

The  Raymond  quarry  is  a  producer  of  white  granite  dimension-stone  products, 
which  are  excavated  from  an  exposed  face  of  an  extensive  mesozoic  granitic 
formation.  The  rock  is  an  even  textured,  fine  grained  moscovlte-blotite 
granite,  grayish  white  in  color,  and  contains  an  extremely  small  amount 


"Sierra  White"  is  the  producer's  trade  name  vor  the  rock  from  this 
quarry. 

**See  Mariposa  Sheet  of  Olaf  P.  Jenkins'  Edition  of  Geologic  Map  of 

California,  Compilation  by  Rudolph  G.  Strand,  U.S.  Geological  Survey,  1967. 
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of  pyrites.  The  chemical  composition  of  the  rock  and  some  of  its  physical 
properties  are  as  follows: 


Chemical  Composition 

Sil  ica  (Si  02) .  73.92% 

Alumina  (Al^O^) . 14.90% 

Ferric  Oxide  ^20^)  and  Ferrous  Oxide  (FeO)  .  .  .  1.95% 

Titania  (T f 02 ) . 0.14% 

Manganese  (MnO) . 0.24% 

Lime  (CaO) . 3.28% 

Magnesia  (MgO)  .  0.68% 

Combined  Water  (^0) . 0.32% 

Soda  (Na20)  and  Potash  (KjO) . 4.51% 

100.00% 

Physical  Properties 

Specific  Gravity  .  2.64 

Bulk  Density  (ASTM  C97-47) . 164.8  lb. /ft. 3 

Absorption  (ASTM  C97-47)  .  0.27% 

Compressive  Strength,  Dry  (ASTM  C 1 70-50)  .  34,800  ps i 

Compressive  Strength,  Wet  (ASTM  C 1 70-5 0 ) .  26,000  ps i 

Modulus  of  Rupture,  Dry  (ASTM  C99-52) .  2050  ps  i 

In  their  excavation  of  the  quarry,  the  operators  have  refrained  from  the 
use  of  explosives  whenever  possible,  and  therefore,  the  rock  mass  is  almost 
entirely  fracture-free. 

While  this  site  is  ideal  geologically,  its  use  presents  other  difficulties, 
primarily  because  of  its  rather  remote  location.  Utilities  and  telephone 
communications  are  available,  but  the  closest  living  accommodations  and 
medical  facilities  are  in  Madera.  The  support  equipment  and  labor  available 
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from  the  quarry  operator  is  rather  limited.  He  will  provide  a  part-time 
electrician  to  provide  power  to  the  test  face,  but  all  other  labor  must 
be  provided  by  West Inghouse .  A  small  machine  shop  and  some  welding 
equipment  are  available  for  Westlnghouse  use,  but  again,  personnel  must 
be  provided  by  Westlnghouse. 
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3.2  SITE  SUPPORT  EQUIPMENT 
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Site  support  equipment  Includes  the  motor-generator  sets  for  the  EBG  power 
supply,  radiation  detection  and  monitoring  instrumentation,  and  radiation 
protection  equipment. 

3.2.1  Electron  Beam  Gun  Power  Suppl les 

Two  portable  engine-generator  sets  have  been  obtained  for  use  in  field 
test  operations.  A  100-kW,  208-volt,  400-Hz,  3-wire  generator  will  be 
used  for  gun  power  and  a  30-kW,  220-volt,  4-wire  unit  will  provide  power 
for  the  auxiliary  equipment. 

3.2.2  Rad lat lcn  Studies 

Wherever  electrons  are  stopped  by  a  target,  x-rays  are  produced.  This  is 

the  mechanism  responsible  for  the  production  of  x-rays  used  for  medical 

and  radiographic  use  and  is  aiso  the  source  of  radiation  in  television 

tubes.  The  stopping  of  electrons  produces  an  x-ray  flux  with  a  continuous 

energy  spectrum,  starting  at  zero  and  reaching  a  maximum  energy  equal  to 

that  of  the  incoming  electrons.  (The  short  wavelength  limit  is  given  by 

^  min  =  12. 4/Eq, where  Eq  Is  the  electron  energy  In  kilovolts.  For  the 

36-kW  electron  beam  gun,  Eq  is  150  kV).  Most  of  the  x-ray  flux  has  an 

eneriy  of  approximately  two-thirds  of  this  value,  or  iOO  kV.  In  addition, 

characteristic  x-rays  of  the  target  material  are  produced.  This  is  mono- 

energetic  "line  radiation".  The  characteristic  K-aipha  line  for  tungsten 

has  a  wavelength  of  X  =  0.209  AU  (corresponding  energy  is  60  kV).  For 

w 

uranium  it  would  be  =  0.126  AU  (99  kV) ;  for  s i 1 i con  »  A  s  j  =  7.111 AU 
(1.75  kV).  The  composition  of  the  target  material  can  be  determined  by 
analysis  of  the  produced  x-ray  spectrum. 

Because  of  the  complexity  of  the  radiation  situation,  one  cannot  give  a 
quantitative  number  for  the  x-ray  flux  produced  and  another  for  the  amount 
of  shielding  required.  To  calculate  the  required  shielding,  one  must 
assume  a  worst  case  and  then  through  the  use  of  empirically  derived  formulas, 
shielding  requirements  can  be  predicted. 
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It  should  be  noted  that  the  assmed  worst  case  results  in  a  conservative 
shielding  design,  since  the  assumed  radiation  energy  level  exceeds  the 
actual  level  by  a  significant  amount.  A  certain  amount  of  "self  shielding" 
also  occurs  as  the  beam  penetrates  into  the  target  material,  but  this  "self 
sh  ielding"  was  neglected  in  the  prediction  of  shielding  requirements. 

3.2.2. 1  Shield  Design.  The  radiation  intensity  produced  by  a  150  kV,  36  kW 
beam  impinging  at  a  point  of  a  flat  surface  of  tungsten  (atomic  no.  7*0. 
according  to  standard  formulas  is  76,500  m  Roentgens  per  hour.  The  intensity 
from  the  point  source  falls  off  with  the  square  of  the  distance  (therefore 
the  m  (square  meter)  term  in  the  above  units).  The  intensity  Is  lower  for 
targets  of  lower  atomic  number,  in  direct  proportion  to  this  number.  For 
silicon  (atomic  no.  14),  which  is  a  major  constituent  of  the  target  rock, 
the  x-ray  intensity  is  therefore  14,500  m  Roentgens  per  hour. 

If  the  operator  is  located  2  meters  from  the  source,  the  unshielded 
radiation  flux  would  be  1/4  of  the  above  number,  namely,  3600  R/hour. 

This  is  the  flux  for  which  shielding  will  be  designed.  A  lead  shield 
gives  the  best  shielding  effect  per  unit  weight  for  the  x-ray  energy 
range  in  which  we  are  working. 

A  concept  design  for  field  test  site  radiation  shielding  has  been  developed, 
based  on  these  intensity  calculations.  The  radiation  shield  design  consists 
of  a  plywood  shed  16  feet  wide  by  about  24  feet  long  by  8  feet  high.  The 
walls  are  lined  with  1/4-inch  lead  in  the  areas  closest  to  the  radiation 
source  and  with  3/16-inch  thick  lead  in  areas  further  removed  from  the 
source.  The  rear  wall  of  the  shield  structure  Is  a  free-standing  baffle, 
which  may  be  removed  with  a  small  crane  to  permit  removal  of  the  gun  and 
support  carriage  without  disassembly  of  the  structure.  Three  leaded  glass 
windows  are  provided  for  observers  and  operators.  The  roof  of  the  shed  Is 
partially  covered  with  heavy  structural  timber  and  a  thin  lead  sheet  which 
offer  protection  from  falling  rock  and  also  reduce  the  intensity  of  air 
scatter  or  "sky  shine"  radiation  to  permissible  levels.  An  artist's  concept 
of  the  shield  structure  appears  in  Figure  3”'. 
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Artist's  Concept 
Structure 


3. 2. 2. 2  Monitoring  and  Detect  ion  Equ ipment .  During  these  initial  field 
tests, extens ive  radiation  intensity  measurements  will  be  made:  (a)  to 
accumulate  more  data  about  the  amount  of  x-rays  actually  produced  with  a 
particular  rock  target,  and  (b)  to  give  absolute  assurance  of  the  operators' 
safety.  Ionization  chamber  type  of  radiation  meters  (calibrated  in 
"Roentgens  per  unit  time"  for  ali  wavelength  of  x-rays)  will  be  used,  as 
well  as  pocket  dosimeters  and  film  badges.  The  gun  control  is  also  equipped 
with  an  audible/visual  Gamma  alarm,  which  is  interlocked  to  the  beam  in 
the  event  of  high  exposure  rates.  (Radiation  safety  codes  require  the  use 
of  film  badges  in  any  radiation  area  where  the  dose  rate  may  exceed  2 
mi i i i roentgens  per  hour.) 
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4.  THEORETICAL  STUOIES 

4.1  PIERCING  WITH  ELECTRON  BEAM 

The  piercing  process  was  fully  described  In  the  First  Quarterly  Report 
(Ref.  3.  P*  >-2).  The  following  discussion  extends  the  previous  work, 
except  that  a  different  scheme  is  used  for  the  computations  of  the  tempera¬ 
ture  distribution  in  the  rock.  The  new  calculations  are  based  on  a  model 
which  will  better  approximate  the  fact  that,  at  the  base  of  the  piercing 
cavity,  the  melt-front  moves  much  faster  than  any  "heat  wave." 

4.2  DEVELOPMENT  OF  THERMAL  STRESSES  WITH  TIME 

Experience  gained  during  the  previous  quarter  concerning  the  rate  at  which 
electron  beams  penetrate  rock  has  been  applied  to  the  calculation  of  the 
resulting  thermal  stresses.  These  stresses,  in  turn,  cause  the  rock  to 
crack.  To  allow  for  eventual  correlation  of  the  results  with  experimental 
data  from  a  9-kW  gun,  the  calculations  assume  an  EBG  that  delivers  9  kW  at 
150  kV  and  Is  1/2  Inch  from  the  rock  face.  The  size  of  the  rock  Is  typical 
of  a  laboratory  sample.  It  is  a  cylinder,  10  Inches  In  radius  and  10  Inches 

high.  The  beam  pierces  the  center  of  one  of  the  circular  ends.  This  gives 

convenient  boundary  conditions  for  the  existing  computer  programs. 

The  calculations  are  me^nt  to  clarify  the  mechanisms  by  which  the  electron 
beam  causes  rock  to  fall.  The  type  of  rock  used  for  the  calculation  Is 
fictitious.  The  cavity  growth  under  an  electron  beam  is  assumed  to  be 
similar  to  that  of  Oklahoma  gabbro  (Ref.  3.  p.  5-2);  the  elastic  parameters, 

similar  to  those  of  Dresser  basalt  (Ref.  4);  and  the  melting  temperature, 

that  of  silica,  l.e.,  1700°  C  (Ref.  5). 

To  provide  the  desired  information  on  the  time  development  of  the  rhermal 
stresses,  six  separate  stress  calculations  were  made;  all  conditions  were 
the  same  except  that  each  successive  calculation  was  made  for  a  longer 
piercing  time.  The  six  elapsed  times  were  5,  10,  20,  40,  80,  and  300 
seconds . 
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4.3  CAVITY  MODEL 

As  mentioned  in  the  First  Quarterly  Report,  the  simple  approximation 
chosen  as  a  model  for  the  typical  cavity  consists  of  a  cylindrical  hole 
with  a  hemispherical  base.  The  cavity  depth  d  is  given  as  a  function  of 
time  t  by  the  equation  (Ref.  3,  P-  4-6), 

f^wo  4^  i  /4 

d  =  12  TT  aH*  1  +  z°  J  ■  zo  ( 1 ) 

where  WQ i s  the  electron  beam  power;  H*,  the  heat  of  "removal"  (partly 
melting,  partly  vaporization,  see  Figure  5“3  of  Ref.  1)  of  the  rock; 
ot,  a  parameter  characterizing  the  spread  of  the  beam  due  to  multiple 
scattering;  and  zQ  ,  the  distance  from  the  original  rock  face  to  the  point 
inside  the  electron  gun  where  the  beam  begins  to  scatter.  Table  4-1  gives 
the  values  of  these  parameters  chosen  for  the  present  calculation  for  the 
9-kW,  1 50— kV  gun. 


Table  4-1 

PARAMETER  VALUES  FOR  DEPTH  VS.  TIME  FUNCTION 
(9-kW,  150-kV  ELECTRON  BEAM  GUN) 


Pa rameter 


Value 


W 

wo 

H* 

a 

z. 


9.0  kW 
223  kJ/in3 
2.78x1 0~3i n"1 
0.941  in 


The  radius  of  the  hole  is  taken  as  its  depth  divided  by  8.4,  a  value 
determined  experimentally  to  be  typical  of  several  types  of  rock  (Ref.  3» 
P.  5-12). 
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4.4  TEMPERATURE  DISTRIBUTION 

To  avoid  lengthy  and  costly  computations  required  to  solve  the  time- 
dependent  heat  flow  equation  with  complicated  boundary  conditions,  the 
following  simplifying  assumptions  are  made  for  temperature-distribution 
cal cul ati ons : 


1.  The  radius  of  the  cylindrical  cavity  remains  constant  while  its 
length  grows. 

2.  The  thermal  diffusivity  is  constant  with  ,..,perature  and  inde¬ 
pendent  of  direction.  The  value  chosen  us  typical  is 

8x10  ^  i n2/sec  (Ref.  6). 


3.  The  radial  temperature  distribution  in  any  one  plane  perpen¬ 
dicular  to  the  z-axis  can  be  calculated  using  the  solution  to 
the  time-dependent  heat  equation  appropriate  for  an  infinitely- 
long,  hot  cylinder  of  constant  temperature,  Tm,  imbedded  in  an 
infinite  medium  at  an  initial  temperature  of  zero  (Ref.  3,  p.  4-9). 
The  solution  is  a  function  of  the  radius  of  the  cylinder,  the 
radial  coordinate  of  the  point  of  evaluation,  the  elapsed  time, 
and  the  thermal  diffusivity. 


4.  The  elapsed  time  used  in  evaluating  the  temperature  distribution 
at  any  plane  z  in  the  rock  depends  on  the  depth  z  of  the  plane, 
since  heating  time  is  longest  at  the  surface  of  the  rock  (z=0) 
and  drops  to  zero  at  the  depth  of  the  hole  (z»d)  and  beyond. 
Specifically,  the  heating  time  is  the  difference  between  the 
piercing  time  needed  to  reach  the  full  depth,  at  which  the  beam 
is  shut  off,  and  the  piercing  time  needed  to  reach  the  depth  z 
of  the  point  of  evaluation.  Using  equation  ( 1 )  yields 


theati  ng^Z  ^ 


2  tt  aH* 
kvt 
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Figure  4-1  shows  this  time  durinj  which  heat  flow  takes  place  in 
any  plane  perpendicular  to  the  z-axis  at  depth  z,  as  though  a 
"snap-shot"  picture  were  taken  at  the  moment  when  the  cavity 
attains  a  depth  of  4  inches.  In  this  and  the  following  figures, 
the  cavity  "depth"  z  is  going  upwards,  because  the  computer 
printouts  were  oriented  and  legible  in  this  position. 


5.  The  value  of  the  cylinder  radius  used  in  the  temperature  function 
at  depth  z  is  taken  to  be  the  radius  of  the  circle  produced  when 
the  plane  defined  by  z  cuts  the  cavity.  In  most  cases,  therefore, 
the  appropriate  radius  is  simply  the  radius  of  the  cylindrical 
portion  of  the  cavity.  However,  at  values  of  z  approaching  the 
depth  of  the  cavity,  the  radius  becomes  that  of  the  circle  formed 
by  passing  the  plane  through  the  spherical  bottom  of  the  cavity, 
and  this  radius  was  used  in  the  calculations. 


Figures  4.2(a)  through  4.2(f)  show  isotherms  plotted  using  the  above  model. 

The  cavity  boundary  is  assumed  to  remain  at  1700°  C;  the  isotherms  are 

o 

shown  in  increments  of  100  C.  Each  figure  represents  one  of  the  six  pierc¬ 
ing  times  considered  in  the  study.  Notice  that  the  isotherms  converge  to 
a  point  at  the  bottom  of  the  hole.  In  this  way,  the  model  approximates 
the  very  steep  temperature  gradient  at  that  point,  which  results  from  the 
advancing  motion  of  the  melt  front.  This  appears  to  be  a  better  approxima¬ 
tion  of  the  real  temperature  distribution  than  the  previously  used  model 
shown  in  Figure  4-2  of  Ref.  i.  The  cavity  depth  increases  much  faster 
than  the  heat  can  diffuse  into  the  rock.  The  extent  of  this  effect  can 
be  estimated  by  assuming  for  the  moment  that  the  speed  of  penetration  is 
constant.  With  this  simplification,  an  upper  limit  can  be  calculated  for 
the  small  distance  that  heat  will  penetrate  into  the  rock  at  the  bottom  of 
the  hole.  The  temperature  rise  at  a  distance  below  the  bottom  of  the 
hole  will  be  given  by  (see  Ref.  7), 


T 


T 


m 


(2) 
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Figure  4-1.  Time  for  which  Heat  Flow  Takes  Place  in  Any  Plane  Perpendicular  to 
Z-axis 
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Figure  4-2.  Calculated  Isotherms  Using  Simple  Model 
with  Cavity  Boundary  at  I700°C 
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where  T^  is  the  melting  temperature  of  the  rock  (the  temperature  of  the 
cavity  wall);  V,  the  speed  that  the  melt  front  is  advancing  (assumed  here 
to  be  constant);  and  D,  the  thermal  diffusivity.  For  V,  the  formula 
in  Ref.  3»  p.  4-6  is  used  to  give 


V  = 


H*  2  tt  a  (d+z0  )3 


(3) 


which  gives  the  penetration  speed  of  the  electron  beam  for  a  given  depth  d. 
(The  parameters  in  equation  (3)  are  the  same  as  those  of  equation  (1)). 
Combining  equations  (1),  (2),  and  (3)  gives 


Az 


AW 

H-'-2  tt  a  r  o 
WQ  \2  TT  olH* 


(4) 


Figure  4-3  shows  plots  of  equation  (4)  giving  the  distance  of  penetration 

of  temperature  into  rock  below  the  end  of  the  cavity  as  a  function  of 

piercing  time  for  several  values  of  temperature  under  the  conditions  listed 

-4  2  o 

in  Table  4-1  wi  th  D  =  8x10  in  /sec  and  T  «*  1700  C.  The  arrows  in 

m 

Figure  4-3  designate  the  times  chosen  for  the  stress  calculations.  As  one 
would  expect,  penetration  of  the  heat  into  the  solid  rock  increases  with 
time  and  for  lower  temperatures.  Contrary  to  the  simplifying  assumption 
made  in  the  above  calculation,  the  speed  V  is  not  a  constant  but  is,  in 
fact,  a  decreasing  function  of  time.  Accordingly,  the  actual  heat  pene¬ 
tration  rates  at  the  base  of  the  cavity  will  be  even  slower  than  those 
given  by  equation  (4). 
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4.5  THERMAL  STRESS  CALCULATIONS 

Thermal  stresses  in  the  rock  were  calculated  using  a  package  of  finite- 
element  computer  programs  whose  details  are  proprietary.  Another  program* 
TEMPMAT,  was  used  to  prepare  the  additional  input  data  temperature  field, 
etc.  for  the  stress  calculations. 

Meshes  for  two  of  the  runs,  No.  2.1  (5  seconds  -  the  shortest  time)  and 
io.  2.6  (300  seconds  -  the  longest  time)  are  shown  in  Figure  4-4.  Note  that 
the  nodal  points  are  spaced  closely  together  in  the  general  direction  of 
the  temperature  gradient  and  are  kept  rather  close  together  in  the  region 
of  greatest  interest,  i.e.,  just  outside  the  heated  region  where  the 
highest  tensile  stresses  are  expected.  Instead  of  starting  at  the  walls 
of  the  cavity,  the  mesh  begins  on  the  700°  C  isotherm,  which  is  taken  as 
the  temperature  at  which  the  rock  becomes  plastic. 

The  program  TEMPMAT  calculates  the  temperature  for  each  mesh  nodal  point 
and  punches  out  the  values,  which  are  then  used  as  inputs  for  the  stress 
calculation  programs.  In  addition,  to  approximate  the  variations  with 
temperature  of  the  elastic  and  thermal  properties  of  the  rock,  TEMPMAT 
makes  use  of  a  material  type-specification  in  the  stress  calculation 
programs  and  assigns  to  each  element  a  value  for  each  of  the  elastic  and 
thermal  parameters  according  to  the  average  temperature  of  the  element 
(see  Table  4-2). 

Figures  4-5  through  4-8  show  the  results  of  the  stress  calculations.  Four 
different  stresses  are  shown:  (1)  azimuthal  tension  (©  stress),  (2)  maxi¬ 
mum  (tensile)  principal  stress  in  the  r-z  plane,  (3)  azimuthal  compression, 
and  (4)  minimum  (compressive)  principal  stress  in  the  r-z  plane. 
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Figure  4-4.  Finite  Element  Triangular  Meshes: 

(a)  Computer  Run  2.1  (5  seconds) 

(b)  Computer  Run  2.6  (300  seconds) 
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Figure  4-6. 


Calculated  Maximum  Principal  Stress 
(Tension)  In  R-Z  Plane  for  Piercing 
Times  of  5,  10,  20,  40,  80,  and 
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Table  4-2 

ELASTIC  AND  THERMAL  PARAMETERS  ASSIGNED  TO  EACH  ELEMENT 


ACCORDING  TO  ITS 

TEMPERATURE 

Temperature 

(°  c) 

Coeff i ci ent  of 
Therm.!  Expansion 

Young's 

Modul us  (psi ) 

Poi sson 
Ratio 

O-iOO 

2.6xl0-6 

145x1 05 

0.2/* 

100-225 

6.0 

141 

0.23 

225-350 

8.2 

135 

0.21 

350-475 

10.2 

126 

0.19 

475-600 

11.2 

107 

0.145 

600-620 

11.7 

83 

0.105 

620-640 

11.8 

64 

0.09 

640-660 

11.9 

46 

0.07 

660-680 

12.0 

28 

0.05 

680-700 

12.0 

9 

0.02 

Note:  These  data  are  similar  to  that  for  Dresser  basalt  In  Table  4-4, 
page  18  of  Ref.  4. 


In  each  of  the  four  cases  there  are  six  plots  representing  the  six  pierc¬ 
ing  times  5,  10,  20,  40,  80,  and  300  seconds.  For  the  plots  of  tensile 
stress,  the  contours  represent  500  to  10,000  psi  In  500  psl  Increments. 

For  the  azimuthal  compressive  stress,  the  contours  represent  0  (dashed  line) 
to  -80,000  ps?  in  -5000  psi  increments.  And,  for  the  compressive  principal 
stress,  the  contours  represent  0  (dashed  line)  to  -76,000  psi  In  -4000  psi 
increments.  The  directions  of  the  principal  stresses  in  the  r-z  plane  are 
more  or  less  parallel  to  the  isostress  contours. 
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4.6  INTERPRETATION  OF  STRESS  FIELDS 

There  are  several  interesting  qualitative  features  of  the  stress  plots  in 
Figure  4-5  through  4-8. 

1.  Compressive  stresses  and  stress  gradients  are  many  times  iarger 
than  tensile  stresses  and  stress  gradients,  i.e.,  typically  by  one 
order  of  magnitude.  (However,  there  is  an  even  wider  disparity 
between  compressive  strength  and  tensile  strength;  hence,  the 
predominant  failure  mechanism  is  failure  under  tension.) 

2.  In  the  neighborhood  of  the  cavity,  the  stresses  are  compressive 
wherever  rock  is  heated.  This  compression  is  caused  by  the 
presence  of  the  surrounding  cold  region,  which  restrains  the 
thermal  expansion  of  the  hot  region. 

3.  Outside  the  heated  region,  the  rock  is  under  tension. 

a.  The  expansion  of  the  hot  region  along  the  axis  of  the  cavity 
tends  to  stretch  the  cold  region  In  the  z  direction.  This 
stress  will  tend  to  cause  cracks  in  a  plane  perpendicular  to 
the  axis  of  the  cavity. 

b.  The  radial  growth  of  the  heated  region  causes  tension  in  the 
azimuthal  direction.  This  stress  will  tend  to  cause  cracks 
in  planes  that  are  parallel  to  and  contain  the  axis  of  the 
cavi ty. 

4.  Azimuthal  tension  is  greater  than  tension  in  the  r-z  plane  at 
depths  down  to  those  comparable  to  the  depth  of  the  cavity.  At 
depths  similar  to  that  of  the  cavity  and  below,  tension  in  the 
r-z  plane  begins  to  exceed  azimuthal  tension.  Such  stresses  can 
produce  failure  configurations  similar  to  those  observed  in  lab¬ 
oratory  experiments,  namely,  a  combination  of  a  vertical  crack 
(caused  by  the  azimuthal  tentlon)  in  a  plane  containing  the  cavity 
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axis  together  with  a  horizontal  crack  (caused  by  the  r-z 
tension),  which  cuts  through  the  bottom  of  the  cavity. 

Figure  4-9  shows  the  best  example  available  of  this  theory.  The 
cavity  in  this  block  of  pink  Jasper  quartzide  is  only  2  inches 
deep  (one  half  the  height  of  the  block),  whereas  in  other  cases 
the  cavity  nearly  penetrated  blocks  of  this  size  before  splitting. 
In  these  latter  cases,  the  horizontal  cracks,  clearly  visible  In 
Figure  4-9,  could  not  develop.  Moreover,  this  quartzite  Is  very 
fine  grain,  homogenous,  and  isotropic.  Note  also  the  fine  cracks 
and  the  lifting  of  the  material  near  the  top  edge  of  the  cavity. 

5.  The  region  over  which  the  stresses  act  expands  with  time,  as  one 
would  expect.  The  depth  of  the  stress  field  exceeds  the  depth  of 
the  cavity,  particularly  for  longer  piercing  times. 

6.  For  each  piercing  time,  the  calculations  show  nontrivial  values 
in  the  cold  region  where  the  tensile  stress  Is  many  times  greater 
than  a  tYBlsal  i^.  strength  of  rosk  (HO3  psl). 
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NOT  REPRODUCIBLE 


Figure  4-9.  Pink  Jasper  Quartzite  Cracked  after  9  Seconds 
of  Piercing  with  9-kW  Beam  (from  Ref.  2) 
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4.7  ROCK  FAILURE 


Although  the  rock  is  expected  to  fail  under  tension,  prediction  of  the 
time,  location,  and  configuration  of  the  fracture  surface  is  quite  compli¬ 
cated.  There  is,  In  fact,  no  unique  tensile  strength  for  a  given  type  of 
rock.  For  example,  Jaeger  and  Cook  (Ref.  8,  p.  140)  make  the  statement: 
"The  tensile  strength  of  rock  Is  more  variable  and  more  Influenced  by 
specimen  size  than  any  other  mechanical  property  of  rock." 

A  useful  interpretation  Is  that  the  effect  of  size  on  tensile  strength 
is  statistical  in  nature.  The  larger  the  volume  under  tension,  the  higher 
the  probability  of  Initiating  a  failure  somewhere  in  the  volume.  This  is 
a  version  of  the  so-called  "weakest  link"  theory.  Jaeger  and  Cook  (Ref.  8, 
p.  186)  discuss  such  a  formula. .on  by  Weibull  (Ref.  9)  which  gives  the 
probability  of  failure  as 


P  ■  '  '  “Xp  [S'n(l-P0)dv]  <5) 

where  dv  i s  an  element  of  volume  and  PQ  Is  the  probability  that  a  unit  volume 
will  fall.  PQ  Is  a  function  of  stress  which,  In  turn,  is  a  function  of  posi¬ 
tion.  For  simplicity,  Weibull  assumes  a  power  law  for  PQ  such  that  P  be¬ 
comes  (in  the  notation  of  Hudson  (Ref.  10)): 

p  -  i  -  «xp  [  -  )m  dv] 

where  m  and  sQ  are  constants  depending  on  the  type  of  rock,  and  s,  a  func¬ 
tion  of  position,  is  the  maximum  principal  tensile  stress. 

Although  Hudson  found  sQ  and  m  to  vary  somewhat  with  volume,  It  is 
instructive  to  use  some  typical  values  that  he  obtained  from  measurements 
on  red  granite,  l.e.,  so^  2000,  and  m»6.  Using  these  values  together 
with  the  calculated  stresses  (Section  4.5),  the  results  shown  In  Figure  4-10 
are  obtained. 
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The  contours  In  Figure  4-10  are  of  the  parameter  I  (r,  s)  defined  by 

r(r-  *)  =(23§o)  2*r 

The  probability  that  a  crack  will  occur  In  a  ring  of  small  cross  section, 

AAsA.rAz  is  then 


P  ■  1  -  exp[  -I(r,  s)Aa] 

The  contours  shown  are  for  I  -  I  to  I  -  10,000;  these  are  derived  from  the 
stress  plot  for  20  seconds  piercing  time  with  a  9  kW,  150  kV  beam. 

From  Figure  4-10,  there  are  apparently  large  regions  that  are  almost  certain 
to  fall  after  only  20  seconds  of  piercing  time.  For  example,  the  co¬ 
efficient  I  reaches  values  of  1000  and  more  near  the  top  of  the  cavity  and 
also  near  the  bottom.  Such  very  localized  regions  where  cracks  are  most 
likely  to  occur  may  trigger  the  all-important  crack  Initiation.  This  may 
happen  after  a  very  short  piercing  time,  because  of  the  steep  temperature 
gradients  produced  by  the  electron-beam  process. 

However,  under  what  conditions  the  cracks  would  then  propagate  and  develop 
is  open  to  question.  For  example,  If  such  localized  stress  begins  to 
develop  early  In  the  piercing  process,  a  small  crack  could  be  produced  and 
relieve  the  stress.  Perhaps  the  neighborhood  surrounding  the  small  crack 
would  subsequently  be  overtaken  by  the  advancing  front  of  molten  rock. 
Further  study  of  this  subject  will  be  necessary  to  determine  the  extent  to 
which  such  a  complicated  sequence  of  events  does  occur.  One  would  also 
like  to  know  to  what  extent  such  a  process  would  alter  the  temperature 
field,  and  at  what  point  the  cracks  will  cease  to  be  localized  and  propa¬ 
gate  far  from  the  cavl ty. 

On  the  basis  of  piercing  tests  and  observations,  one  would  ~xpect  that  the 
larger  volumes  heated  to  lower  temperatures  and  farther  from  the  melt-front, 
as  well  as  farther  from  the  top  surface,  are  responsible  for  propagating  the 
large  fractures,  which  do  occur. 
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5.  LABORATORY  TESTS 

No  laboratory  tests  were  performed  during  the  second  quarter. 
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